Recent calculations and analyses of O star spectra 
have revealed discrepancies between theory and observa- 
tions, and between different theoretical calculations, for 
the strength of optical Hei singlet transitions. We inves- 
tigate the source of these discrepancies. Using a non-LTE 
radiative transfer code we have undertaken detailed test 
calculations for a range of O star properties. Our princi- 
pal test model has parameters similar to those of the 09V 
star, 10 Lac. We show that the discrepancies arise from 
uncertainties in the radiation held in the Is 2 1 S-ls2p 1 P° 
transition near 584A. The radiation field at 584A is in- 
fluenced by model assumptions, such as the treatment of 
line-blanketing and the adopted turbulent velocity, and by 
, the Feiv atomic data. It is shown that two Feiv transi- 
'w? ' tions near 584A can have a substantial influence on the 
, strength of the He I singlet transitions. Because of the dif- 
ficulty of modeling the He I singlet lines, particularly in 
stars with solar metalicity, the He I triplet lines should be 
preferred in spectral analyses. These lines are much less 
sensitive to model assumptions. 



00 



Key words, stars: atmospheres - stars: early-type - line: 
formation - radiative transfer 



(N 

in 
o 

O 

S3 

Pl,: 

6 



X 



Astronomy Sz Astrophysics manuscript no. 4489ms 


February 5, 2008 


(DOI: will be inserted by hand later) 





On the sensitivity of Hei singlet lines to the Feiv model atom 

in O stars 

F. Najarro 1 , D. J. Hillier 2 , J. Puis 3 , T. Lanz 4 , and F. Martins 5 

1 Instituto de Estructura de la Materia, Consejo Superior de Investigaciones Cientficas, CSIC, Serrano 121, 28006 
Madrid, Spain 

2 Dept. of Physics and Astronomy, University of Pittsburgh, 3941 O'Hara St., Pittsburgh, PA 15260, USA 

3 Universitats-Sternwarte Miinchen, Scheinerstrasse 1, 81679 Munich, Germany 

4 Department of Astronomy, University of Maryland, College Park, MD 20742, USA 

5 Max Planck Institiit fur Extraterrestrische Physik, Postfach 1312, 85741 Garching, Germany 

Received DD MMMM 200Y/ Accepted DD MMMM 200Y 
Abstract. 



1. Introduction 

With the advent of non-LTE line-blanke ted m odel atmo- 
sphere calculations (Hubeny and Lanz Il995t Hill ier k . 
Miller Il998t Pau ldrach et al. l200lt Grafener et al. l2002t 
Puis et al. 120051) considerable advances in the analysis of 
massive stars have been made. There has, for example, 
been a significant downward revisio n in the effe ctive tem- 
perat ure of O stars (Marti ns et al. I2002L 120051: Crowther 
et al. l2002t Herr ero et al. l2002t Bian dii fc Garcia l2002t 
Garc i a fc Bi anchi 12004 Repolust et al. 120041: Massey et al. 
I2004L 1200.51) . Moreover, the discrepancy between spec- 
troscopic and evol utiona ry masses (Groenewegen et al. 
Il989t H errero et al. I1992T) ha s been greatly reduced (e.g., 
Herrero l200l Repolust l20ol|) . 

A key diagnostic of the effective temperature in O stars 
is the ratio of He 1 lin e stre ngths relative to tho se of He 11 
(e.g., Kudritzki et al. ll983l Herrero et al. l2000|) . However 
detailed comparisons with observations show discrepan- 
cies larger than expected. Further, theoretical modeling 
reveals sensitivities of some lines to the adopted tech- 
niques and assumptions while there are also discrepancies 
between results obtained using different codes. 

From very early on it was realized that in extended 
atmospheres the ls2p 3 P° and ls2p 1 P° states become 
overpopulated, and that t his affects the streng th of the 
He 1 lines (e.g., Wellmann Il952al Il952bi Il952ct Ghobros 
Il962[) . This is termed the dilution effect , and was more re- 
cently discussed by Voels et al. | fl989) . With the advent 
of non-LTE line-blanketed model atmospheres it was ex- 
pected that the new model atmospheres would be capable 
of modeling the full Hei spectrum. However, as noted by 



Hillier et al. | ((2003J) , it was not possible to get a simulta- 
neous fit to all He 1 lines. Further, in that study of two O 
stars in the SMC, it was noted that the Hei singlet lines 
were more sensitive to model details than the He 1 triplet 
line 4473 A. As additional studies have been undertaken 
it has become clear that while model calculations often 
have difficulty fitting the He 1 singlet transitions they give 
satisf actory fits for th e Hei triplet lines (e. g., Pu is et al. 
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120051 Repo lust et al. 120051 Ma rtins et al. 
et al. l2005l Mokiem et al. 120051: Lanz et al. 

More r ecently , comparisons between CMFGEN (Hillier 
and Miller Il998|) and TLUSTY (Hubeny and Lanz ll995|) . 
and between CMFGEN and FASTWIND (Puis et al. 
2005), have shown inconsistent predictions, in some pa- 
rameter ranges, for the s trength of the Hei singlet 
lines (e.g., Puis et al. l2005|h The disagreement between 
CMFGEN and TLUSTY was somewhat su rprisi ng since 
previous comparisons (e.g., Hillier & Lanz l200ll Bouret 
et al. 120031) had shown excellent agreement. However the 
comparisons were limited to a few models, and it is has 
become apparent that the discrepancies only occur for cer- 
tain parameter ranges. Given the different assumptions 
in the codes, different atomic data, and different tech- 
niques, it is not immediately obvious which code is giv- 
ing the correct answers. Agreement between observations 
and CMFGEN was improved in some cases, for example, 
by adding additional species (e.g., Argon, Neon, Nickel 
and Calcium) , and by reducing the turbulent velocity used 
in the atmospheric structure/population calculations from 
20 to 10 km s" 1 (Martins et al. l2005|) . 

In this paper we discuss the analysis of Hei lines in 
O stars, and the inconsistencies between triplet and sin- 
glet lines, between observations and models, and between 
results obtained using different atmospheric codes. The 
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complicated behavior of the Hei singlet line strengths is 
shown to be a consequence of complicated radiative trans- 
fer effects in the neighborhood of the He I Is 2 1 S - Is 2p 1 P° 
transition at 584A. In particular, it is shown that two Feiv 
lines, which closely coincide in wavelength with this tran- 
sition influence the population of the ls2p 1 P° level, and 
thus the strength of Hei singlet lines that couple to the 
Is 2p 1 P° level. Indeed, a closer examination of the discrep- 
ancies reveals that the lines that are most discordant be- 
tween different model calculations are not the He I singlet 
lines per se, but rather the transitions involving Is 2p 1 P°. 

That line- line interactions, and the precise details of 
line blanketing in the far UV, might affect optical line 
diagnostics is not surprising. In planetary nebula the fluo- 
rescenc e of O ill lines by He n Lya (the Bowen mechanism, 
Bowe n ^934) is well documented (see, e.g., Osterbrock 
Il989h . In WN stars many of the N lines are affected by 
continuum fluoresc ence p rocesses, and hence sensit ive to 
blanketing (Hillier flQSsf) . Further, Schmutz l|l997h sug- 
gested that photon loss from the Hen resonance transi- 
tion, by overlapping metal lines, could significantly affect 
the ionization structure in W-R stars and hence help to 
solve the problem of how to accelerate W-R winds. The ef- 
fect was investigated in the modeling the WC comp onent 
of the WR+O binary 7 Velorum (De M arco l2000lh and 
found to be important. Martins | l|20A4h investigated the 
same mechanism in O star models (i.e., the influence of 
metal lines near He 11 Lya) and also found that the effects 
were significant. 

2. The interaction of He 1 with Fe iv 

A simplified Grotrian diagram for the He 1 atom (singlets 
only) is shown in Fig. For the present discussion the 
most important transition is the ls 21 S- ls2p 1 P° transi- 
tion at 584 A. Key transitions, often used for diagnostic 
purposes, are indicated. Overlapping the 584 A transition 
are several Feiv transitions. The data for the relevant 
Feiv tr ansitio ns used in CMFGEN, taken from Kurucz 
& Bell l|l995h . is provided in Table 1 — only the two 
most important transitions are shown. Wavelengths for 
the Feiv transitions were calculated using the known 
energy levels given in the NIST Atomic Spectra Database 
(http://physics.nist.gov/PhysRefData/ASD/index.html), 
and are believed to be accurate. 1 The velocity separation 
of the two Fe IV lines from the Hei A584 is given, a s is th e 
oscillator strengths computed by Becker & Butler (l995). 
Thes e show significant differences to those of Kurucz & 
Bell (|l995|) . and indicate that the oscillator strengths for 
these transitions must be regarded as uncertain. 

In O stars Hei 584 A is often optically thick, so that 
photons undergo a number of scatterings within the line 

1 The energy values for the relevant lower levels in the NIST 
table are quoted to 0.1 cm -1 , while the upper levels are quoted 
to an accuracy of 0.01 cm -1 . An accuracy of 0.1cm -1 in the 
energy difference between the two levels corresponds to less 
than 0.2 km s -1 , and thus it is unlikely that the wavelengths 
are significantly in error. 
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Fig. 1. Simplified Grotrian diagram for the singlet states 
of He 1. The A584 transition plays a key role in determin- 
ing the population of the ls2p 1 P° level, and hence the 
strength of transitions that involve that level. Although 
not shown, another important line often used as a diagnos- 
tic (A4389; ls2p 1 P°-ls 5d : D) also involves the ls2p 1 P° 
state. 

before escaping. Because the optical depth is large, a high 
population of the ls2p 1 P° state is maintained enhanc- 
ing the strength of absorption lines ending on this state. 
However, as the photons scatter other mechanisms come 
into play. In particular, they can be absorbed by the over- 
lapping Feiv lines. As the upper states of these levels have 
alternative decay routes, the absorbed photons can be re- 
moved from the 584 transition, acting as a drain, and caus- 
ing the Is 2p 1 P° population to be lowered. This leads to a 
weakening of the transitions whose lower state is Is 2p 1 P°. 
Clearly the importance of this effect in model atmosphere 
calculations will depend on whether the explicit interac- 
tion between the Hei resonance transition and the Feiv 
lines is taken into account, the accuracy of the Fe IV wave- 
lengths, oscillator strengths and branching ratios, and on 
other atmospheric parameters such as the adopted micro- 
turbulent velocity. It is also obvious that the importance 
of the effect will depend on the actual stellar parameters 
(i.e., T ef f, log g, M). 

In this study we concentrate specifically on the influ- 
ence of the two Feiv lines on the Hei line strengths. It 
should be stressed, however, that other Feiv lines, and 
lines due to other species (e.g., Cav, Fevi) could also po- 
tentially have an influence (albeit generally weaker) on the 
Hei singlet line strengths. 

3. Results 

In order to quantitatively explore the influence of the Fe iv 
lines on the Hei lines we consider the following model: 
T eff = 33,560K, log 5 = 3.9, L = 7.13 x 1O 4 L , R* = 
7.8 R Q , M = l.OxlO^Mgyr- 1 , = 1400 kms" 1 . The 
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Ion 


Low level 1 


Upper level 


A(A) 


AV^kms" 1 ) 2 


A ut 


f 3 

J lu 


,4 
J lu 


He i 


Is 2s *S 


Is 2p 'P 


20586.9 




1.73 x 10 u 


0.329 




He I 


Is Is 4 S 


Is 2p 1 P° 


584.334 




1.78 x 10 9 


0.274 




Feiv 


3d 5 4 F 9/2 


3d 4 ( 3 G)4p 2 H° /2 
3d 4 ( 3 G)4p 4 H? /2 


584.368 


17 


6.82 x 10 7 


0.00349 


0.00251 


Feiv 


3d 5 2 D3 5/2 


584.397 


32 


4.22 x 10 s 


0.0288 


0.00264 



1. Level designations for Feiv are from the NIST Atomic Spectra Database. The D3 designation is used to denote one of the 
three 3d 5 2 D terms. 

2. Velocity shift of Feiv transition relative to He 1 A 584 tr ansition. 

3. Feiv oscillator strengths are from Kurucz & Bell 



4. Feiv oscillator strengths are from Becker & Butler (1995 



parameters of the model are similar to models we are us - 
ing to investigate the 09V star 10 Lac (Lanz et al. [2006) . 
Since the mass-loss is low, the important parameters for 
the present study are T e f f and log g only. The hydrostatic 
structure was initially obtained from a TLUSTY model 
but has been subsequently modified as the parameters 
were changed. Below 5kms _1 hydrostatic equilibrium is 
achieved to generally better than 5%. Unless otherwise 
stated, we used a fixed Doppler width of 13.5 km s -1 in 
the atmosphere calculations, while a turbulent velocity of 
lOkms -1 was adopted for the profile calculations. The 
N(He)/N(H) ratio is taken as 0.1 while we adopted solar 
abund ances using the compilation of Grevesse & Saval | 
l)l998|) . CNO abundances were reduced sli ghtly in accord 
with the recent results by Asplund et al. 1 l)2005|) . 

3.1. He I singlet sensitivities: f values, abundances, 
and turbulent velocities 

In order to test the sensitivity of the models to the oscilla- 
tor strengths of the two relevant Fe iv transitions we ran 
a detailed CMFGEN model using the tabulated / values, 
and then models in which the / values were reduced by fac- 
tors of 2, 5, and 10. The results for Hei A4923 (vacuum), 
A5017, A20587, and A4473, are shown in Figures[l|3||3|and 
[^respectively. The three singlet lines show significant vari- 
ations. For the transition with ls2p 1 P° as the lower level, 
the line is strongest in the model with the lowest / val- 
ues. As noted earlier Fe iv transitions drain photons from 
the Is 2 4 S - ls2p 1 P° transition, lowering the equilibrium 
ls2p 1 P° population. When the oscillator strength is low- 
ered the drain is weakened, and the Is 2p 1 P° population is 
enhanced. The A4923 transition, which is coupled directly 
to the ls2p 1 P° state, shows a large variation (Fig. [2J, 
while the A5017 line (Fig. which has 2s *S as its lower 
state, shows smaller variations. Another transition, He I 
A6680, which also has ls2p 1 P° as its lower state, shows 
similar or even more sensitive behavior than A4923. The 
near IR transition, Hei A20587 (Is 2s ^-ls 2p 1 P°), also 
shows very significant variations. However, in this case, 
the absorption profiles are strongest when the / values are 
largest. This is to be expected, since for this transition the 
ls2p X P° state is the upper level. The difficulty in model- 
ing the A20587 line in early type stars has be en dis cussed 
extensively, for example, by Najarro et al. I|1994l) . The 
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Fig. 2. Predicted line profiles for Hei A4923 (2p 1 P°-4d 
4 D). The solid curve show the predictions with the stan- 
dard model, while the other curves show the profiles when 
the relevant Feiv / values are reduced by a factor of 2 
(dashed curve), a factor of 5 (dot-dash), and a factor of 
10 (dot). 

problem with A20587 has not been fully solved with mod- 
els accounting for bla nketing. In their analysis of WR147, 
Morris et al. | l)2000|) found that this line was one of the 
very few lines for which a discrepancy still appeared. They 
obtained a much too deep He 1 20587 line, consistent with 
the singlet problem presented here. The He 1 A20587 line 
is very important for analyses of O stars using H and K 
band infrared data (Repolust et al. |2005|) . and hence it is 
important that its formation be understood. 

The triplet line, Hei A4473, is essentially identical in 
the 4 models (Fig. [SJ). Similar statements also apply to 
other Hei triplet lines (e.g., A5877). Limited tests with 
TLUSTY and FASTWIND confirm the changes found us- 
ing CMFGEN. 

Another obvious way to reduce the "drain" is to lower 
the Fe abundance. Indeed, as seen from Fig. El a reduction 
in the Fe abundance by a factor of 2 gives a similar en- 
hancement in line strength in A4923 as does the factor of 
2 reduction in the Fe line strengths. The correspondence 
is not exact since a lower Fe abundance reduces line blan- 
keting and hence causes slight changes in the temperature 
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Fig. 3. As for Fig. but for Hei A5017 (Is 2s ^-ls 3p 
1 P°). The line strengths have changed significantly, but 
the change is less than for transitions ending on Is 2p 1 P°. 




20582 20584 20586 20588 20590 20592 

Fig. 4. As for Fig. but for Hei A20587 (2s x S-2p 1 P°). 
Because the Is 2p 1 P° state is the upper level of this tran- 
sition, the variations in the strength are in the opposite 
sense to the optical lines. 
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Fig. 5. As for Fig. but for Hei A4473 (2p 3 P°-4s 3 D). 
The triplet line strengths are not affected by changes in 
the Feiv line strengths. 
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Fig. 6. Predicted line profiles for Hei A4923 (2p ^Md 
1 D). The solid curve show the predictions with the stan- 
dard model, the dashed curve for a factor of 2 reduction 
in the Fe abundance, and the dot-dash curve for a factor 
of 2 reduction in the relevant Feiv / values. 



structure of the model atmosphere. As before, the A4473 
profiles are very similar in the three models. The sensi- 
tivity of the singlet lines to the Fe abundance explains 
why Hillier & Lanz l)200l|) achieved excellent agreement 
when comparing models for SMC O stars where the Fe 
abundance is only 0.2 solar. 

As a final test of the sensitivity of the Hei line 
strengths we examined the effect of the microturbulent 
velocity used in the atmosphere structure calculations to 
determine the level populations. In Fig.[7|we show the pro- 
files for A4923 computed using the same turbulent velocity 
for the profile calculations (lOkms -1 ) but computed us- 
ing fixed Doppler widths of 10 and 13.5 km s -1 for the 
determination of the atomic populations. As might be ex- 
pected, A4923 is stronger when computed using the lower 



turbulent velocity. With the larger turbulent velocity there 
is better overlap between the Fe iv and the He I line, and 
hence a more effective drain. 

The use of lOkms -1 for population calculations in O 
stars is not unusual. Detailed analysis suggest that this 
turbulent velocity is fairly typical, with possible lower 
values for dwarfs, and somewhat higher values for super- 
giants. The precise physical cause/meaning of microtur- 
bulent velocities derived from profile fitting is unclear. 

We have also performed test calculations for other 
"typical" O star parameters. Similar effects are seen to 
occur, although the size of the changes depends strongly 
on spectral type. The mechanism does not appear to be 
operative for B stars, and in a model of a late O supergiant 
(T eff = 28, 400 K, log g = 3.2, M = 1.65 x lO^Moyr" 1 , 




4921 4922 4923 4924 4925 4926 
A(A) 

Fig. 7. Predicted line profiles for Hei A4923 (2p 1 P°-4d 
1 D). The solid curve shows the predictions with the stan- 
dard model while the dashed curve shows the effect of 
a lower turbulent velocity in the calculation of the level 
populations. 
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A(A) 

Fig. 8. Predicted line profiles for Hei A4923 (2p ^"-Ad 
1 D) in a model for a late O supergiant. The solid curve 
shows the predictions with the standard model, while the 
other curve shows the profile when the relevant Fe iv / 
values are reduced by a factor of 10. 

L = 5.4 xlO 5 ^©) the influence of the Fe IV lines was rather 
small, as illustrated in Fig.|Hl Conversely, the influence in 
an early O dwarf (T e// = 41, 300 K, log 5 = 4.0, M = 
2.0 x 10~ 6 M©yr-\ L = 2.35 x 10 5 i Q ), where the Hei 
lines are weak, is relatively strong (Fig.EJ). Interestingly, 
in a mid-0 supergiant (T e ff — 36,200K, logg = 4.0, 
M = 9.5 x lO- 7 M yr _1 , L = 2.35 x 10 5 L Q ), a reduction 
of a factor of 10 in the two Fe iv oscillator strengths can 
cause He I A4923 to change from being weakly in emission 
to a relatively strong absorption line. 

In principle, / values could be empirically derived from 
a model comparison with observation. However, given the 
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Fig. 9. Predicted line profiles for Hei A4923 (2p 1 P°-4d 
*D) in a model for an early O dwarf. The solid curve shows 
the predictions with the standard model, while the other 
curve shows the profile when the relevant Fe IV / values 
are reduced by a factor of 10. 

problem with turbulent velocities, and the sensitivities to 
the Fe abundance and other stellar parameters, we refrain 
from providing actual values. However, it does appear that 
both th e Kur ucz & Bell (1995) values, and the Becker & 
Butler l|l995|) values, are too large. 

Taken as a whole, these findings explain the differ- 
ent results obtained between TLUSTY, CMFGEN, and 
FASTWIND. While TLUSTY and CMFGEN use similar 
data for Feiv the models are not always computed us- 
ing the same assumptions, such as which species are in- 
cluded and the adopted turbulent velocity. In CMFGEN, a 
smaller turbulent velocity leads to a longer computational 
time, and so a larger turbulent velocity is often adopted. 
In general the choice has very little effect of the popula- 
tion determinations but for a few lines, as shown here, the 
adopted value can be important. The adopted turbulent 
velocity was also shown to directl y affec t the Hen A4687 
line in B supergiants (Evans et al. l2004|) . In FASTWIND 
blanketing is treated using mean opacities and emissiv- 
ities, and consequently the detailed interaction between 
the Is 2 1 S-ls2p 1 P° lines at 584A and the Feiv lines is 
not taken into account, although work is in progress to 
rectify this. Since all three codes give similar answers for 
the He I triplet lines, the triplet lines are clearly preferred 
for spectral analyses. 

As noted by the referee, Herrero et al. | l|2002|) were 
able to obtain a good fit to the He 1 A4473, and A4389 
lines using a model with T e /r = 35, 500 ± 500 K and 
\ogg — 3.95 ±0.1. This can be explained by noting that 
in 10 Lac it appears that the neglect of the influence 
of the two Feiv transitions on the Hei singlet transi- 
tions is a reasonable approximation. Finally, it should 
be noted, that FASTWIND consistently gives tempera- 
tures, for O dwarfs, that are 1500 K higher than CMFGEN 
(and TLUSTY). The reason for the discrepancy is un- 
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known but is being investigated. Mokiem et al. 
also obtained similar fit parameters to that of Herrero 
et al.(T e// = 35,000 ± 850 K and \ogg = 4.03 ± 0.13) 
using numerous FASTW1ND models and a sophisticated 
genetic algorithm. However, they derived a turbulent ve- 
locity of 15.5 ± 4kms _1 which is inco nsisten t with that 
indicated by the metal lines (Lanz et al. l2006[) . again sug- 
gesting that possible systematic effects in the models are 
being neglected (see additional discussion in Lanz et al. 

4. Conclusion 

Earlier work on the analysis of massive O stars revealed in- 
consistencies in the predictions of He I line strengths with 
observations, and between different codes. We have shown 
that these inconsistencies arise because of the different 
treatment of line-blanketing in the neighborhood of the 
Hei ls 21 S - ls2p 1 P° transition at 584A, due to differ- 
ent treatments regarding the interaction of the Hei and 
Feiv model atoms, and due to different adopted model 
Feiv atoms. The adopted microturbulent velocities and 
iron abundance also play a direct role in determining the 
population of the ls2p 1 P° state, and hence on the Hei 
singlet line strengths. 

The accurate prediction of Hei lines, which involve 
ls2p 1 P°, is difficult, and possibly not feasible with the 
current atomic data available for Fe iv. We argue strongly 
that the Hei triplet lines should be weighted the most 
heavily when performing spectral analyses. Inconsistencies 
between model fits of the He I singlet lines and triplet lines, 
and between those singlets connected to Is 2p 1 P°(e.g., 
A4923) and other singlet lines (e.g., A5017; ls2s 1 S- 
ls3p 1 P°) are indicative of significant transfer effects in 
Hei A584. 

Extensive modeling by one of us (Najarro) suggests 
that the Kurucz Fe iv oscillator strengths are too large but 
given the uncertainties in the models, and the influence 
of the turbulent velocity, it is not feasible to determine 
reliable empirical values. 

We have shown how subtle radiative transfer effects 
can significantly influence the strength of the He I singlet 
transitions. It is possible that similar effects may occur for 
other diagnostic lines, especially if the levels involved in 
such transitions are coupled to low lying levels whose pop- 
ulation may be effected by overlap between strong transi- 
tions and metal lines. 

Interestingly, the subtle effect of the Fe iv lines on the 
ls2p 1 P° population may contribute to the dilution ef- 
fect. In early unblanketed O star models, triplet lines were 
consistently too weak compared with observations , espe - 
cially in stars of class II, lb, and la (Voels et al. |l989). 
This was interpreted as a problem with the triplet lines, 
and was believed to arise from the neglect of extension 
effects. Intrig uingly, all the singlet lines analyzed by Voels 
et al. I l)l989|) involve the ls2p 1 P° state. While it is now 
known that blanketing is much more important than was 
previously suspected, the use of spherical models, the al- 



lowance for stellar winds, and the influence of line blan- 
keting (as in FASTWIND) has not alleviated the problem 
with H 1 A4473 in giants and supergian ts of t ype 06 and 
later (Massey et al. l2005l Repolust et al. l2004h . Physically, 
and with the new generation of models, it is much easier 
to understand the inconsistencies in the He I line strengths 
as arising from problems in computing the population of 
the ls2p 1 P° state and hence the singlet line strengths, 
rather than the triplet line strengths. 

This work highlights the importance of performing 
analyses with independent codes, and with performing 
systematic analyses of all available lines. In this way dis- 
crepancies, sensitivities, and biases can be determined and 
investigated. The work also highlights the importance of 
understanding the physics of line formation if one is to 
derive accurate abundances. This is especially true in sit- 
uations where only one or two spectral lines are available 
for analysis — a situation not uncommon in O stars. 
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